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1 Introduction 



Magnetic fields are ubiquitous in the Universe, being observed in nature on scales from 
planetary size to galaxy cluster size |Q, 0- In galaxies and galaxy clusters, the typical 
strength is of order a few //Gauss, which is thought to be produced by dynamo action on 
a seed field. In galaxies the dynamo timescale is roughly a rotation period, 10® Yr, and a 
simple calculation [|l| based on the age of a typical galaxy shows that the seed field must 
have been about 10^^° Gauss, or perhaps less in the currently favoured models with a 
cosmological term [Q. 

There is no shortage of ideas for generating this seed field. The more conventional 
astrophysical explanations are based on a Biermann battery operating at the era of reion- 
isation (see e.g. and the references therein). There are more speculative ideas based 
on various models of inflation [Q, phase transitions and primordial black holes [§. 

All these mechanisms have the common feature of producing stochastic, homogeneous 
and isotropic magnetic and velocity fields which can be characterised by their power 
spectra and characteristic initial scales. Our interest here is to try and make model- 
independent statements about the evolution of the magnetic fields once they are gener- 
ated. This article, which is based on Ref. 0, studies the evolution of a stochastic magnetic 
field generated at a phase transition, such as the confinement transition in QCD at t ~ 1 
sec, or the electroweak symmetry-breaking transition at t ~ 10~^^ sec. It therefore falls 
into the category of decaying 3D MHD turbulence, which has been studied before in the 
MHD community [0, |ll], . Most directly comparable to our work, Biskamp 

and Miiller studied the energy decay in incompressible 3D magnetohydrodynamic tur- 
bulence in numerical simulations at relatively high Reynolds number, and in a companion 
letter studied the scaling properties of the energy power spectrum. 

We focus here on the transfer of magnetic energy from small to large scales, as ne- 
cessitated by the conservation of magnetic helicity. This is important for a primordial 
magnetic field to reach a large enough scale with sufficient amplitude to be relevant for 
seeding the galactic dynamo [jl6|. 

We perform 3D simulations both with and without magnetic helicity, starting from 
statistically homogeneous and isotropic random initial conditions, with power spectra 
suggested by cosmological applications. We find a strong inverse cascade in the helical 
case, with the coherence scale of the field growing as t^'^, and equivocal evidence for a 
weak cascade when only helicity fluctuations are present. In the helical case we also find 
a self-similar power spectrum with an approximately behaviour at high k. We find 
decay laws for the magnetic and kinetic energies of t""'^ and t~^'^, respectively, in the 
helical case, and t^^'^ for both in the non-helical case. 
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2 Evolution of magnetic fields in the early Universe 



A convenient benchmark is to assume that the field is created on the horizon scale with 
a power spectrum /c", taking all the energy in the Universe, and that it is subsequently 
completely frozen into the plasma |]TB|. If one takes the epoch of creation to be either the 



QCD transition or the electroweak transition, one finds that the RMS fiuctuations on the 
scale of a protogalaxy are roughly 

B{to, L = O.SMpc) < I jg-20-4.5r!^Q 

Thus we see that for causal fields (n > 2) there needs to be some amplification on large 
scales for there to be any interesting seed field for the galactic dynamo. 

In fact, the field is not strictly frozen into the plasma, and in order to calculate 
observable effects, we must determine how the field scale length ^ and magnetic energy 
Em evolves. There are various scaling arguments which have been put forward: for 
example, for ideal MHD (infinite conductivity) and with no helicity Olesen (and later 
Son |18[, Field and Carroll [0 and Shiromizu [^) argued 



e(t)~t2/("+^), (1) 

where t is conformal time. The effect of having a conserved helicity modifies this scaling 
law to g 0, |1|] 

m ~ t'/\ Em ~ t-'/\ (2) 



Early numerical experiments with a shell model of the full MHD equations ||2T]] suggested 
^ ~ t^''^^, and focussed attention on the possibility of an inverse cascade, in which power 
is transferred locally in /c-space from small to large scales. 



3 MHD equations 

The MHD equations in an expanding Universe are most conveniently expressed in terms 
of conformally rescaled fields B, u and dissipation parameters z/, t], and in the gauge 
A^ = r]V-u 11,1]. 

The matter and radiation in the early Universe is modelled as an isothermal com- 
pressible gas with a magnetic field, which is governed by the momentum equation, the 
continuity equation, and the induction equation, written here in the form 

-u- Vu-c^Vlnp+:?^ + ^fv^u+^VV-uV (3) 

p p \ S / 

-u ■ V In p - V ■ u, (4) 

u X B + r^V^A, (5) 



du 

'dt 
dlnp 

dt 
dA 

'dt 
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where B = V x A is the magnetic field in terms of the magnetic vector potential A, u is 
the velocity, J is the current density, p is the density, /x is the dynamical viscosity, and 77 
is the magnetic diffusivity. 

In the ideal limit p = rj = 0, there is a conserved quantity in addition to the energy, 
which is the magnetic helicity Hu, given by 

Hm= [ A B d^x. (6) 



Helicity is known to be important in dynamo theory [^, |2^, where turbulence is driven. 
We shall also be able to confirm its importance in decaying turbulence. 

An important dynamical quantity is the magnetic Reynolds number Reu = Lv/rj, 
where L and v are the typical length scale and velocity of the system under consideration, 
because it measures the importance of the non-linear term in the equation for the magnetic 
field. In the early Universe, Reu can be very large. This is often taken to mean that 
the magnetic field is frozen into the plasma, and the scale length of the field increases 
only with the expansion of the Universe. This is in general untrue, because turbulence 



can transfer energy to different length scales pT[. Not only can there be the usual direct 



cascade of energy from large to small scales, but also an inverse cascade, increasing the 



overall comoving correlation length [23 . 



4 Plasma properties in the early Universe 

We are considering the time before recombination, when the plasma in the early Universe 
has a relativistic equation of state p = p/3, and therefore a sound speed = l/y/S. The 
number species contributing to the pressure gradually decreases with the temperature T, 
until T ~ rrie, the electron mass, when only photons and neutrinos remain relativistic. At 
that point the number of charge carriers in the plasma reduces by a very large factor: for 
T ^ rrie, the electron number density rie is approximately equal to the photon number 
density n^, while for T <C rrie n^/n^ = ns/n^ — 10~^°, where Ub is the baryon number 
density. 

The transport properties of the plasma are determined from the mean free path /^fp of 
the relevant particles involved in the transport of the quantity of interest, which typically 
is Imfp ~ l/a| loga|T |2^, where a is the fine structure constant. From this we can infer 
that the conductivity is 

{T/a\ loga| T ^ rrie 

{ne/nr^){ml/T)/a T <^ 711^ , (7) 

x{nB/n^){T/meyTe-^ T < T^,, 

where Tdoc is the temperature at photon decoupling, roughly 1 eV. 
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The viscosity parameter v is given by 



_ / / l/a^l logajr T 3> me , . 

^ = ^/P ~ I (n^/n,)K/r2)/ar T « m, 

We can estimate the importance of non-hnear terms relative to diffusion from the 
Reynolds numbers of the plasma, the hydrodynamical Re = vLjv and magnetic Rej\,/ = 
vLa. Upper bounds are obtained by assumed a fluid flowing at the speed of light on a 
scale somewhere between the mean free path and the particle horizon ct: 

(a, Xjc? \oi a) ^ (Re, ReM) ^ 10^^(T/GeV)-^ (9) 



5 3D MHD simulations of decaying turbulence 



We solve these equations numerically with a code |^ using a variable third order Runge- 
Kutta timestep and sixth order explicit centred derivatives in space. All our runs are 
performed on a 120^ grid, and we use periodic boundary conditions, which means that 
the average plasma density (po) = Po is conserved during runs. Here po is the value of the 
initially uniform density, and the brackets denote volume average. 

We use natural units c = 1, and set the unit of length by setting measure fci = 1, 
where k\ is the smallest wave number in the simulation box. Hence the box has a size of 
-^^Box = 27r. The scale factor is fixed by setting po = 1? and B is measured in units of 
^/popoc, where po is the magnetic permeability. We define the mean kinematic viscosity 
V = p/po- The sound speed Cg = as appropriate for a relativistic fluid. 

The equations are not quite those for a relativistic gas in the early universe ||21|| . 
However, we have checked that our results change little when using the true relativistic 
equations in the low velocity limit. 

The initial conditions appropriate for the early Universe are to take both the velocity 
and magnetic fields to be homogeneous and isotropic Gaussian random fields drawn from 
a power-law distribution with a high wavenumber cut-off, determined by the physical 
mechanism generating the power spectrum. The power spectra, defined by Pm(^) = 
(BJ^ ■ Bk), and P\{k) = (uj^ ■ u^) is taken to have the form 

Puik) = AMk^exp{-{k/k,)'), Py{k) = Ayk"' exp{-{k/k,Y). (10) 

note that in the plots it is the shell-integrated energy spectra, Em,v = 4:7ik^ x |-Pm,v(^), 
which are shown. 

Note also that causality demands that n > 2 and m > In the simulations 



presented we took the power laws to be the lowest consistent with causality. We also 
chose kc = 30, unless specified otherwise. The initial magnetic energy was taken equal to 
the kinetic energy, and had the value 5 x 10~^ in all runs. 



4 



We are also able to control the helicity in the initial conditions, ranging from identically 
zero to maximal, where maximal helicity is defined as saturating the inequality 

\HM{k)\<2k-^EM{k). (11) 

It should be noted that if one chooses the initial field with no constraint, there are helicity 
fiuctuations present, around a mean of zero. It is possible to completely eliminate the 
fiuctuations but in practice they are quickly regenerated by the evolution of the field. 



6 Results 

In all runs the mean kinematic viscosity v and the resistivity r] were chosen to be equal 
with values between u = rj = 5 x 10"'^ — 5 x 10~^. 



6.1 Magnetic energy spectrum 

In Fig. m we show magnetic energy spectra Eyi{k) for runs with unconstrained and maxi- 
mal initial magnetic helicity. Turning first to maximal initial helicity, one sees that after 




Figure 1: Magnetic energy spectrum Euik) for a run with unconstrained (left) and maxi- 
mal (right) magnetic helicity. u = rj = 5 x 10~^. The times shown are 0, 1.0, 4.6, 10.0, 21.5 
and 46.3. At low wavenumbers k the energy spectrum Euik) increases with time. 



a short initial direct cascade, there is energy transfer towards lower fc, and the peak of 
the power spectrum moves to the left, signifying a growth in the coherence scale of the 
field. This has been called an inverse cascade, although it should be noted that the term 
"cascade" carries the implication of a local interaction in fc-space, which is not necessarily 
true p9|]. In fact, in simulations where the initial conditions contain power in only a single 



wavenumber, the k"^ spectrum is very quickly established at low k 
the power transfer to large scales is non-local in fc-space [p6 



30], indicating that 
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For wavenumbers above the peak, one sees a decrease in power with quite a differ- 
ent form than is present in the initial conditions, which can be fitted approximately by 

When only helicity fluctuations are present in the initial conditions, the peak in the 
power spectrum still moves to low k, but there is very little transfer of power. It is perhaps 
surprising to see any transfer at all, as there are arguments which connect the evolution 
of the coherence length with helicity conservation p, |1^ . In an attempt to test this idea 
we also did a run with identical parameters but zero initial helicity, but saw essentially 
identical behaviour, with a small increase in power at small scales. 

In our runs with ii^K = -^^m initially, the kinetic energy spectrum shows no evidence 
of an inverse cascade at any scale. However, when the initial velocity distribution is zero 
the kinetic spectrum grows on all scales initially and in the low wave number region the 
energy continue to grow even after the high wave number modes start to decay. 

6.2 Coherence length evolution 

One length scale is the magnetic Taylor microscale Lt, defined as -Brms/-^rms- In Fig- 0we 
show Lt against time t for a run with maximal initial helicity. The asymptotic behaviour 



10° 10^ 10^ 

t 

Figure 2: Time evolution of the magnetic Taylor microscale for with maximal initial 
magnetic helicity. u = t] = 5 x 10~^. The straight line indicates the power law oc 

of the length scale is seen to grow approximately as Lt ~ t^'^- 

In runs with non-helical initial conditions the growth of the magnetic Taylor microscale 
is slower: we find approximately Lt ~ t"'^. 
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6.3 Energy decay laws 



Figure ^ shows the magnetic energy Eu{t) and the kinetic energy EY^{t) with maximal 
initial helicity. It is seen that the asymptotic decay rate for Eyiif) is approximately 




Figure 3: Time evolution of the magnetic energy Eyi{t) and the kinetic energy E^{t) in 
the case where there is initial magnetic helicity. u = rj = 5 x 10~^. The straight lines 
indicate the power laws oc and oc t^^'^ respectively. 

The kinetic energy also decays with a power law behaviour at late times: -Ek(^) ~ t^^'^- 
In runs without initial helicity the decay rates of £'m(^) and -Ek(^) are approximately 
the same, close to t^^-^. 

6.4 Magnetic Reynolds number 

The Reynolds numbers in our simulations are evaluated using the magnetic Taylor mi- 
croscale. In most our simulations we typically obtain Reynolds numbers of the order of 
100 — 200. In Fig. ^ we show the magnetic Reynolds number of the run with maximal 
helicity whose power spectrum is shown in Fig. |I|. Note that for the second half of the 
run the magnetic Reynolds number is approximately constant. 

6.5 Self-similarity in magnetic power spectrum 

We make the following ansatz for the energy spectrum 

EM{k,t) = m-''gMm, (12) 

where ^ is the characteristic length scale of the magnetic field, taken to be the magnetic 
Taylor microscale defined above, and g is a parameter whose value is some real number. 
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Figure 4: Magnetic Reynolds number for the run of Fig. || 
Figure ^ shows ^{t)'^Euik,t) versus the scaled variable k^(t), with q = 0.7, plotted for 
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Figure 5: The magnetic scaling function guik^) described in the text, equation (|T2D, 
against kC,, with power laws (kC,)'^ and (kC,)'"^'^ for comparison. 

different values of time t. The data is seen to collapse onto a single curve given guikC,): 
thus the magnetic field evolves in a self-similar manner. 



7 Discussion and conclusions 

The most directly comparable simulations of decaying 3D MHD turbulence were carried 



out by Biskamp and Miiller They found similar results, the magnetic field evolved 

self-similarly, with a power-law behaviour at high k. However, their power law was k~^^^, 
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much less steep than our A;"^'^. There were a number of differences between their and 
our simulations: they were able to achieve larger Reynolds numbers, both by having 
larger grids, and by using hyper-diffusivity (a magnetic diffusivity term). However, 
we believe that the real difference is due to their initial cut-off scale being significantly 
larger, at kc = 4. We have performed a run with larger initial length scale, kc = 5. In this 
case the magnetic energy spectrum develop into an approximate k~^^^ law at late times. 
However, this occurs only after the peak of the spectrum has left the simulation box. 

There have also been renormalisation group (RG) analyses looking for an inverse 
cascade in driven MHD turbulence [2C, 28 1. In particular Berera and Hochberg p8 



saw 



no evidence for an inverse cascade. However, it is not clear that the results are directly 
comparable, firstly because we are considering freely decaying turbulence, and secondly 
because RG analysis can only give information about late times when the system is in 
equilibrium with the driving force. What we have referred to as an inverse cascade is. 



even in the driven case [Q, a time-dependent phenomenon characterised by a bump in 
the power spectrum travelling to smaller wave number. 

In conclusion, we find good evidence from our numerical simulations that helical 
stochastic magnetic fields show an inverse cascade (in the sense explained above), and 
that even if only small helicity fluctuations are present initially, there is still weak inverse 
cascade. 

We have determined growth laws for the magnetic and kinetic energies Em and Ek- 
In the helical case, Em ~ t~'^''^ and Ek ~ t~^'^, which means that that there is no 
equipartition of energy. This is because the extra constraint of helicity conservation forces 
the magnetic held to transfer power to larger scales rather than allow it to be dissipated. 
The importance of helicity is borne out by the fact that in the non-helical case, we find 
^ ~ and Em Ek ^ t'^'^ ■ 

Length scales in the magnetic field increase as t^'^ in the helical case, but slightly slower 
in the non-helical case, t^'^. Note that these growth laws disagree with all theoretical 
predictions to date [T^, |l8l , which give t^^^ in the helical case and t^/^ for our power 



spectrum in the non-helical case. 

Helical magnetic fields are found to evolve in a self-similar way, with a scaling function 
gM{z) ~ at large k, where p = 2.5 for Re ~ 10^. Note that this is significantly 
different from both the Iroshnikov-Kraichnan and Kolmogorov spectra, and k~^^^, 

respectively. 

A good theoretical understanding of these scaling laws is required before the evolu- 
tion of magnetic fields in the early Universe is properly understood, as a small error in 
the exponent makes a large error in the prediction of the magnetic field strength when 
propagated over many orders of magnitude in time. For example, Vachaspati's contri- 
bution to these proceedings |^ assumes a growth law of t^^^ in the length scale, based 
on a simple argument invoking helicity conservation [Q, |T8|, |l9l, to obtain seed fields of 
an interesting strength from the electroweak transition, which is quite different from our 
observed growth law of 



9 



ACKNOWLEDGEMENTS 



This work was conducted on the Cray T3E and SGI Origin platforms using COSMOS 
Consortium facihties, funded by HEFCE, PPARC and SGI. We also acknowledge com- 
puting support from the Sussex High Performance Computing Initiative. MH thanks 
NORDITA for hospitality while this work was completed. 



References 

[1] Ya.B. Zeldovich, A. A. Ruzmaikin and D.D. Sokoloff, Magnetic Fields in Astrophysics 
(Gordon & Breach, New York, 1983); A. A. Ruzmaikin, A.A.Shukurov and D.D. 
Sokoloff, Magnetic Fields in Galaxies (Kluwer, Dordrecht, 1988). 

[2] P.P. Kronberg, Rep. Prog. Phys. 57, 325 (1994). 

[3] D. Biskamp, Nonlinear Magnetohydrodynamics (Cambridge Univ. Press, Cambridge, 
1993). 



[4] A. C. Davis, M. Lilley and O. Tornkvist, Phys. Rev. D 60 (1999) 021301 ||arXiv:astro- 
ph/990402l . 



[5] N.Y. Gnedin, A. Ferrara and E.G. Zweibel Astrophys.J. 539 (2000) 505-516 
||arXiv:astro-ph/000l0661 



[6] M. S. Turner and L. M. Widrow, Phys. Rev. D 37 (1988) 2743.; B. Ratra, 
Ap. J. 391, LI (1992); A. Dolgov, Phys. Rev. D 48 (1993) 2499 [ |arXiv:hep 



ph/9301280|| .; M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. Lett. 



75 (1995) 3796 ||arXiv:hep-th/9504083|| .; D. Lemoine and M. Lemoine, Phys. Rev. 



D 52 (1995) 1955.; O. Bertolami and D. F. Mota, Phys. Lett. B 455 (1999) 96 
arXiv:gr-qc/9811087| |; K. Dimopoulos, T. Prokopec, O. Tornkvist and A. C. Davis, 



||arXiv:astro-ph/0108093l 



[7] C.J. Hogan, Phys. Rev. Lett. 51, 1488 (1983); J.M. Quashnock, A. Loeb and D.N. 
Spergel, Ap. J 344, L49 (1989); B. 1. Cheng and A. V. Olinto, Phys. Rev. D 50 
(1994) 2421; G. Sigl, A. V. Ohnto and K. Jedamzik, Phys. Rev. D 55 (1997) 4582 
||arXiv:astro-ph/961020l| ; M. M. Forbes and A. R. Zhitnitsky, Phys. Rev. Lett. 85 



(2000) 5268 ||arXiv:hep-ph/0004051 



[8] A. Vilenkin and D. Leahy, Ap. J 248, 13 (1981); Ap. J 254, 77 (1982). 

[9] M. Christensson, M. Hindmarsh and A. Brandenburg, Phys. Rev. E 64, 056405 
(2001) | ^rXiv:astro-ph/00ll32l |. 



[10] M. Hossain, P. Gray, D. Pontius, W. Matthaeus and S. Oughton, Phys. Fluids 7, 
2886 (1995). 



10 



[11] H. Politano, A. Pouquet and P.L. Sulem in Small-Scale Structures in Fluids and 
MHD, edited by M. Meneguzzi, A. Pouquet, and P.L. Sulem, Springer- Verlag Lecture 
Notes in Physics Vol. 462 (Springer- Verlag, Berlin, 1995), p. 281. 

S. Galtier, H. Politano and A. Pouquet, Phys. Rev. Lett. 79, 2807 (1997). 

M. Mac Low, R.S. Klessen, A. Burkert, and M.D. Smith, Phys. Rev. Lett. 80, 2754 
(1998). 

D. Biskamp and W.C. Miiller, Phys. Rev. Lett. 83, 2195 (1999). 
W.C. Miiller and D. Biskamp, Phys. Rev. Lett. 84, 475 (2000). 
M. Hindmarsh and A. Everett, Phys. Rev. D58, 103505 (1998). 
P. Olesen, Phys. Lett. B398 321 (1997) ||arXiv: astro-ph / 9708004|| . 
D. T. Son, Phys. Rev. D 59 (1999) 063008 [|arXiv:hep-ph/9803412 |. 
G. B. Field and S. M. Carroll, Phys. Rev. D 62 (2000) 103008 ||arXiv:astro- 



ph/9811206 



T. Shiromizu, Phys. Lett. B 443 (1998) 127 jarXiv: astro-ph / 98 103391 . 
A. Brandenburg, K. Enqvist, and P. Olesen, Phys. Rev. D54, 1291 (1996). 
K. Subramanian and J. D. Barrow, Phys. Rev. D 58 (1998) 083502 ||arXiv:astro- 



ph/9712083| . 



A. Pouquet, U. Frisch, and J. Leorat, J. Fluid. Mech. 77, 321 (1976). 

M. Meneguzzi, U. Frisch, and A. Pouquet, Phys. Rev. Lett. 47, 1060 (1981). 

J. Ahonen and K. Enqvist, Phys. Lett. B 382 (1996) 40 ||arXiv:hep-ph/9602357|| ; 
J. Ahonen, Phys. Rev. D 59 (1999) 023004 ||arXiv:hep-ph/980l434|| ; G. Baym and 



H. Heiselberg, Phys. Rev. D 56 (1997) 5254 ||arXiv:astro-ph/9704214 



A. Brandenburg, ApJ 550, 824 (2001) ||arXiv:astro-ph/000618e 



R. Durrer, T. Kahniashvih and A. Yates, Phys. Rev. D 58, 123004 (1998) 
||arXiv:astro-ph/9807089 . 



A. Berera and D. Hochberg, |arXiv:cond-mat/0103447 . 
MH would like to thank Jens Niemeyer for clarification of this point. 
M. Christensson, M. Hindmarsh and A. Brandenburg, in preparation. 
T. Vachaspati | arXiv: astro-ph /0 1 1 1 1 24 



11 



